Structure and dynamics of topological defects in a glassy liquid on a negatively curved 

manifold 



Frangois Sausset, 1 '^ Gilles Tar j us, 2 '[j] and David R. Nelson 3 '^ 

'institut de Physique Theorique, CEA, CNRS URA 2306, F-91191 Gif sur Yvette, France 
2 LPTMC, CNRS-UMR7600, Universite Pierre et Mane Curie, 
boite 121, 4 Place Jussieu, 75252 Paris cedex 05, France 
3 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 

(Dated: March 17, 2010) 

We study the low-temperature regime of an atomic liquid on the hyperbolic plane by means of 
molecular dynamics simulation and we compare the results to a continuum theory of defects in 
a negatively curved hexagonal background. In agreement with the theory and previous results on 
positively curved (spherical) surfaces, we find that the atomic configurations consist of isolated defect 
structures, dubbed "grain boundary scars", that form around an irreducible density of curvature- 
induced disclinations in an otherwise hexagonal background. We investigate the structure and the 
dynamics of these grain boundary scars. 



I. INTRODUCTION 

More than a century after Thompson's considera- 
tions about the organization of repulsive charged par- 
ticles in spherical shells [1] , the problem of ordering and 
minimum-energy configurations of particles on curved 
surfaces has attracted renewed interest A theory 

has recently been proposed which derives ground-state 
particle configurations on 2-dimensional (2D) curved 
manifolds from the effective free-energy of relevant in- 
teracting topological defects, disclinations and disloca- 
tions [3l [6]-[9] . It has been mainly applied to "spherical 
crystals", i.e. defected triangular tilings on the surface 
of a sphere, for which it predicts the appearance of stable 
"grain boundary scars" formed by strings of dislocations 
emanating from an irreducible curvature-induced discli- 
nation and terminating in the crystal. Such defect struc- 
tures are energically forbidden in flat (Euclidean) space 
and result from the nonzero curvature of space. This pre- 
diction, as well as more specific ones concerning either the 
structure or the dynamics of the defects, have been con- 
firmed in both simulations [5j and experiments |10H12] 
on spherical crystals. 

The purpose of the present work is to test the theoreti- 
cal predictions on 2D manifolds of constant negative cur- 
vature. Negative curvature brings in new features com- 
pared to the more familiar positive curvature of spherical 
geometry. First, a homogeneous space of constant neg- 
ative curvature can be of infinite extent: in 2D, this is 
the hyperbolic plane H 2 . As a result, a thermodynamic 
limit and bona fide phase transitions may be envisaged, 
contrary to what occurs on a sphere. Secondly, the hy- 
perbolic plane as a whole cannot be embedded in 3D Eu- 
clidean space so that physical realizations of negatively 
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curved manifolds in soft and hard condensed matter do 
not have constant curvature. They can be negatively 
curved locally, in the form of a saddle point, or in a more 
extended way, as on the inner wall of a torus [5] . They 
can also be of arbitrary spatial extent with a negative 
Gaussian curvature only on average, such as the "mini- 
mal surfaces" found for example in bicontinuous phases 
of amphiphilic bilayers |13j . 

Here, we consider an equilibrium system of atoms 
interacting with a circularly symmetric pair potential 
on the hyperbolic plane, which we investigate through 
molecular dynamics simulation. This model has been 
introduced and already studied in the frustration-based 
approach of glasses [U HH [TS] , of supercooled liquids and 
the glass transition [T6Tll8| . and of jamming [19j[20]. In 
the present context of low-energy configurations of par- 
ticles and associated topological defects, we focus on the 
model at low temperature, in a regime where the spatial 
extension of the local hexagonal/hexatic order charac- 
teristic of the liquid has saturated due to the frustrating 
effect of space curvature [IS]. We find that the atomic 
configurations consist of rare defect structures in an oth- 
erwise 6-fold coordinated background. At the tempera- 
tures at which equilibration is still possible on the simu- 
lation time scale, the system is disordered at long length 
scales and appears as a very slowly relaxing "glassy" liq- 
uid. Interestingly, the defect structures are organized 
much like what has been found in spherical geometry 
and predicted by the continuum theory of defects EJ- 
9 . They form grain boundary scars emanating from an 
irreducible number of disclinations forced in the medium 
by the topology of the embedding manifold. By chang- 
ing the magnitude of the curvature, we also confirm, at 
least in a semi-quantitative way, several predictions of 
the theory. 

The rest of the article is organized as follows. In section 
II, we introduce the atomistic model and we shortly de- 
scribe basic notions and specific features coming with the 
hyperbolic geometry as well as the way we handle them 
in the computer simulations. In section III, we review 
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the continuum theory of defects and apply it to nega- 
tively curved 2D manifolds. Section IV is devoted to the 
simulation results and their comparison with theoretical 
predictions. Finally, we give some concluding remarks in 
section V. 



II. MODEL, SIMULATION, AND BASIC 
NOTIONS 



We consider a one-component atomic liquid embedded 
in the hyperbolic plane H 2 . The latter is a homoge- 
neous surface of constant negative curvature K = —n 2 , 
the associated hyperbolic metric being obtained from 



ds 2 



' sinh dip with (r, ip) the appropriate 



polar coordinates [5TJ [55] . The atoms interact through 
the Lennard- Jones potential 



v(r) = 4e 



r 



(1) 



where r is the geodesic distance between atom centers. 
We have investigated this model by means of molecular 
dynamics simulation. The simulations have been per- 
formed at fixed density pa 2 ~ 0.85 in the NVE ensem- 
ble, and several curvature magnitudes, parametrized by 
kg, have been considered. We have restricted our study 
to "small" curvatures (k<t < 0.2) for which the local ar- 
rangement of the atoms is the same as in flat (Euclidean) 
space, i.e. hexagonal. Due to the curvature-induced frus- 
tration, hexagonal tiling of the whole plane is however 
forbidden and an irreducible density of defects must be 
present. For larger curvatures (kct > 1), local heptago- 
nal, octagonal, etc., atomic arrangements become pre- 
ferred in the liquid and, for specific ranges of values, this 
may lead to unfrustrated (quasi) long-range ordering in 
regular lattices with coordination numbers z > 7 at low 
temperature |15j . 

Hyperbolic geometry requires extending the usual 
methods and algorithms used in Euclidean space in order 
to account for boundary effects (which are never neg- 
ligible compared to bulk contributions), parallel trans- 
port of vectors (which is not globally defined on curved 
manifolds), and visualization (since H 2 cannot be em- 
bedded in 3D Euclidean space). Details on the proce- 
dures are given elsewhere [TTJ [TS] . We just mention here 
that, to study bulk behavior, we have implemented pe- 
riodic boundary conditions [53]- Quite generally, peri- 
odic boundary conditions can be envisaged from both a 
"geometrical" and a "topological" perspective. From a 
geometrical point of view, a periodic boundary condi- 
tion amounts to replicating a primitive cell containing 
the system of interest in order to generate a tiling of the 
whole space. From a topological one, a periodic bound- 
ary condition implies providing a rule to pair the edges 
of the primitive cell, which leads, after identification of 
the paired edges, to a "quotient space" that can be rep- 
resented as a multi-connected compact manifold. The 
primitive cell together with the edge pairing rule define 



a "fundamental polygon" [53]. On the hyperbolic plane, 
there is an infinite number of possibilities. A given fun- 
damental polygon is then characterized by an area, or 
equivalently by the genus of the associated compact man- 
ifold (see below Eq. Q), a number of sides, which for 
regular polygons and for a given genus fixes the asso- 
ciated tiling, and a pairing of its sides. In the present 
study, we have used two different boundary conditions. 
The first one corresponds to an octagonal primitive cell 
whose quotient space can be topologically described as a 
compact manifold of genus g — 2 (a "2-hole torus"); the 
associated tiling of the hyperbolic plane H 2 is denoted 
by {8,8}, meaning that 8 octagons meet at each ver- 
tex (more generally, the Schlafli notation {p, q} denotes 
a tiling in which q regular polygons with p edges meet at 
each vertex). The other periodic boundary condition is 
based on a a polygonal cell with 14 edges and a quotient 
space with genus g = 3 ( "a 3-hole torus" ) ; the associated 
tiling of H 2 is a {14, 7} one. Choosing these two different 
conditions allows us to check that the observed behavior 
of the liquid does correspond to bulk behavior with neg- 
ligible finite-size effect: the area of the 14-gon is indeed 
twice bigger than the area of the octagon whereas, as will 
be shown below, the measured quantities are essentially 
identical in both cases. 

Visualization, on the other hand, is obtained through 
the Poincare disk model [211 E2]- This is a projection 
of the whole hyperbolic plane H 2 onto a unit disk. The 
projection is conformal (angles are preserved) but not iso- 
metric (distances are deformed and shrink as one reaches 
the disk boundary). Others representations of H 2 exist, 
each having its own advantages: the Poincare half-plane, 
the Klein model, the hyperboloid of revolution endowed 
with the Lorentz metric (ds 2 = dx 2 + dy 2 -dz 2 ) [5T1I55] . 
An often used way to represent H 2 in real space is the 
so-called pseudosphere: it does not need any projection, 
has a constant negative curvature, and is of infinite ex- 
tent. Its shape is defined by the surface of revolution of 
a tractrix [22] around its asymptote and can be roughly 
seen as a trumpet horn of infinite height. We have how- 
ever chosen the Poincare disk model because it is well 
suited for the visualization of particles and defects. 

In 2D space, topological defects are point-like. A mi- 
croscopic definition of the elementary defects, which are 
disclinations in a reference hexagonal crystal, is possible 
through a Delaunay triangulation of the atomic config- 
urations. The procedure can be generalized to the hy- 
perbolic plane, for which it has been shown to be unique 
when the ensemble of points formed by the atomic cen- 
ters is dense enough [53], a condition always satisfied in 
our system. From the Delaunay construction, one can 
then determine the coordination number of each atom in 
the liquid. When na < 0.2 and at low temperature T, 
most atoms are 6-fold coordinated, which corresponds to 
a local hexagonal arrangement, and defects are located 
on atoms whose coordination number differs from 6: for 
instance, a 7-fold coordinated atom corresponds to a neg- 
ative disclination of topological charge — 7r/3 whereas a 
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5-fold one corresponds to a positive disclination of charge 
7r/3; disclinations of higher charges ±qn/3 with an inte- 
ger q > 2 are very rare at low T . For illustration, two 
configurations at high T (above the ordering transition 
T* in flat space) and low T (significantly below T*) are 
shown in Fig. [I] In this work, we only focus on defect 
structures at low temperature, when the spatial exten- 
sion of the local hexagonal order (which grows as T de- 
creases) has saturated as a result of curvature-induced 
frustration [18]. 

In 2D space, there are two important identities, the 
Euler-Poincare relation and the Gauss-Bonnet theo- 
rem [5TJ [55] . The former only involves topology and for 
N points on a compact surface that is triangulated via 
the Delaunay construction, it reads 



N 



(6-c)=x = -2(fl-l) 



(2) 



where c is the average coordination number of the points 
and where we have used the definition of the Euler char- 
acteristic of the surface in terms of the genus, x = 
— 2(g — 1). On the other hand, for a surface £ of total 
area A embedded in the hyperbolic plane H 2 of curvature 
K = —k 2 , the Gauss-Bonnet theorem states that 



f d 2 x^g~(^)K = 2K X , 



(3) 



where <?(x) is the determinant of the metric tensor at 
point x, which in polar coordinates is simply equal to 
[sinh(Kr)/ft] 2 . The curvature being constant, Eq. (3| can 
be rewritten as 



A k = 47r(.g - 1) 



(4) 



which combined with Eq. ^ gives the deviation from 6 of 
the average coordination number in terms of the density 
of particles p = N/A and the curvature: 
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(5) 



One can apply Eqs. ^ and ([5 1 to the elementary cell of 
the tiling of H 2 generated by the periodic boundary con- 
ditions (this cell corresponds to the fundamental poly- 
gon defining the periodic boundary condition, see sec- 
tion |n|. Eq. ([2| shows that, neglecting disclinations of 
charge magnitude higher than tt/3 |27j . i.e. particles 
with coordination number c < 4 and c > 8, the cell must 
contain an excess of exactly 12(g — 1) negative 7-fold 
disclinations (the remaining 7-fold disclinations exactly 
compensate the number of 5-fold disclinations). Thus, an 
excess of 12 such disclinations is present for the octagonal 
cell {g = 2) and of 24 for the 14-gonal cell (g — 3). This is 
to be contrasted with the result on a sphere where an ex- 
cess of 12 positive 5-fold disclinations is required. In ad- 
dition, Eq. Q shows that the area of the elementary cell 
depends on both the curvature and the type of bound- 
ary condition (through the genus g). Note then that the 



excess density of 7-fold disclinations is independent of 
the specific fundamental polygon chosen for the periodic 
boundary condition and only depends on the curvature. 
It is indeed equal to 12(g - 1)k 2 /[4w(g - 1)] = 3k 2 /ir. 



III. CONTINUUM THEORY OF DEFECTS IN 
HYPERBOLIC GEOMETRY 

In this section, we summarize the continuum theory 
developed in Refs. [5] [BHS] • We more specifically apply it 
to negatively curved manifolds, which is appropriate for 
our atomic model on H 2 with periodic boundary condi- 
tions. The formalism directly deals with defect degrees of 
freedom and starts with the elementary ones, the discli- 
nations, from which all other topological defects can sub- 
sequently be built. Six-coordinated particles are treated 
using continuum elastic theory [3]. At low temperature, 
the elastic free energy of an arbitrary disclination density 
on a 2D manifold S of constant negative curvature — k 2 
is given by [3 EH] 



F =~ J J d 2 x^) d 2 y^g(j) x 



( S (x) 



JLl 
A2 l*y 



(6) 



>(y) 



where Y is the Young modulus of the hexagonal crys- 
tal in flat space, 1/A 2 is the inverse Laplacian-squared 
operator, and s(x) the disclination density, defined as 



S (x) 



(7) 



where qi € Z* (i.e., all integers except zero) is the re- 
duced topological charge of the ith disclination located 
at Xj (in practice at low T, we only need to consider 
qi = ±1) and Nd is the total number of disclinations 
in S; finally, E core ^ is a disclination core (free) energy, 
renormalized by thermal fluctuations 3J. Note that in 
the present problem, Eq. Q is the free energy for the 
elementary cell of the tiling of H 2 associated with the 
periodic boundary condition. Provided interactions be- 
tween cells are weak, the total free energy of the whole 
system is just this quantity multiplied by the (macro- 
scopic) number of cells. This is of course different from 
spherical geometry where the whole system is finite. The 
disclination density must therefore satisfy a "topological" 
generalization of the electroneutrality constraint, namely, 



d 2 xy/glxj (s(x) + k 2 ) =0, 



(8) 



even at finite temperature. By using Eq. ([7| and the 
Gauss-Bonnet theorem, Eq. Q, the above expression 
leads to a sum rule on the total topological charge, 
namely 



N d 



'12(5-1). 



(9) 



(a) 



(b) 



FIG. 1: (Color online) Atomic configurations in the Poincare disk representation of H 2 . Two temperatures above 
and below the ordering temperature T* ~ 0.75 in flat space are shown: (a) T/T* ~ 2.4, (b) T/T* ~ 0.52. From the 

Delaunay triangulation, atoms are colored according to their coordination: black for 6 (hexagonal), red for 7 
(negative disclination of charge — 7r/3) and blue for 5 (positive disclination of charge +7r/3); other colors correspond 
to disclinations of higher charge magnitude. Here, k<t = 0.05 and the octagonal periodic boundary condition (pbc) is 

used. 



Not surprisingly, one then obtains the same result as pre- 
viously derived on the basis of the Euler-Poincare rela- 
tion: there must be an excess of 12(g — 1) disclinations 
with elementary charge q = — 1 (per elementary cell) at 
low T. 

The presence of the 12(g — 1) localized irreducible de- 
fects only imperfectly screens the curvature which is uni- 
form throughout space. The resulting elastic strain en- 
ergy could then potentially be reduced by introducing 
dislocations that represent tightly bound pairs of +1 and 
— 1 disclinations. The amplitude of the stress generated 
by one of the irreducible, unpaired, —1 (i.e. 7-fold) discli- 
nations is controlled by an effective disclination charge 
at a geodesic distance r of the disclination center. This 
quantity is given by 



eff(r) 



2tt 



dr' Sinh(Kr,) 



(10) 



[7- 6cosh(«x)] 



and vanishes at the critical radius r c = 
k- 1 cosh -1 (7/6) ~ 0.57k" 1 . 

It was shown in Ref. [51 [BJ for spherical geometries 
that there is a range of dimcnsionless curvatures no~ for 
which the strain energy associated with the stress field 
generated by the irreducible disclinations is minimized by 
the presence of strings of dislocations that radiate from 
each of the irreducible disclinations. This can easily be 



transposed to hyperbolic geometry. These strings of dis- 
locations, whose Burgers vector is perpendicular to the 
string axis, terminate in the crystalline background at a 
distance of the order of the critical radius r c . These struc- 
tures have been dubbed "grain boundary scars" [51 151 HP]. 
The typical number of dislocations per irreducible discli- 
nation can be estimated at zero temperature by assum- 
ing that a scar contains m strings aligned along geodesies 
emanating from the disclination center. It has then been 
argued [3] that dislocations along a geodesic are roughly 
separated by a distance t(r) ~ m<j/\s e tf(f)\, so that 
the number of dislocations per irreducible disclination 
is given by 



N d ~m I dri(r) -1 - est 
'o 



-est 



TV 



7 cosh 1 



(w)~ 



(11) 



-est + 0.3998 (kct)~ x , 



where est denotes a constant coming from subdominant 
terms which depend on microscopic details (such as the 
Young modulus and the core energies). In contrast, the 
contribution proportional to (kct) -1 is universal for a 
given geometry and underlying type of crystalline order. 

The above estimate does not take into account possible 
interactions between grain boundary scars. In the same 
approximation, one can derive from Eq. ^ the (free) en- 
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Irreducible 
disclination 




FIG. 2: (Color online) Sketch of a grain boundary scar 
with 2 arms of dislocations formed by 5 — 7 (blue-red) 
dipoles attached to a central 7-fold (red) irreducible 
disclination. The glide motion takes place along the 
Burgers vector as indicated by the arrow. 



ergy of a boundary scar. Upon assuming that the strings 
radiating from an irreducible disclination and forming a 
given scar do not interact, the energy can be written as 
the sum of dislocation-dislocation pair energies, which 
are always repulsive, and of a disclination-dislocation en- 
ergy, which for a dislocation in a geodesically aligned 
string at a distance r from the disclination center is given 
by 



F Dd (r) 



Y 
12 



{ko)- x 



sinh(Kr) In 



cosh(Kr) — 1 
2 



1 + cosh(/tr) 



(12) 



where we have assumed that the magnitude of the Burg- 
ers vector is given by <r, the lattice spacing of the under- 
lying crystal. This expression can be obtained from that 
derived for spherical crystals [9] by replacing by ire. 
Fdcl{t) is attractive up to a distance of order It 
is the competition between the energy gain coming from 
this attraction and the energy cost due to the repulsion 
between dislocations and to the dislocation core energies 
that results in stable grain boundary scars, at least for 
sufficiently small dimensionless curvature ko~. 

We now turn to the dynamics of the defects at low tem- 
peratures. We focus on the structural relaxation time r 
at low but finite T. We again consider a situation in 
which there are 12(g — 1) noninteracting grain bound- 
ary scars. A simple estimate of the structural relaxation 
time is provided by the time scale for the diffusion of the 
most loosely bound dislocation in a scar to another scar. 
Climb motion being energetically prohibitive, the motion 
of a dislocation proceeds via thermally activated glide 
motion (i.e. parallel to the Burgers vector). This takes 
place in a potential energy that combines the Peierls po- 
tential due to the underlying lattice and the interaction 
potential with both the central disclination and the other 
dislocations of the scar [12] . The main contribution to the 
activation energy barrier along the dislocation motion is 
due to the disclination-dislocation interaction. The latter 



increases as the distance between the outermost disloca- 
tion and the central disclination increases along the glide 
direction (which, we recall, is perpendicular to the aver- 
age direction of the string: see Fig. [2]). A crude estimate 
for the barrier is obtained from a harmonic approxima- 
tion to Fud in Eq. ( 12 ) for a transverse displacement of 
the dislocation. For small curvature magnitude k<j and 
after taking into account that the outermost dislocation 
is originally at a distance ~ re -1 of the central disclina- 
tion, this leads to 



F Dd {K-\y) - F Dd ( K -\0) ~ CYnay 2 + 



(13) 



where C > is a constant of order unity and y is the 
distance along the geodesic that is perpendicular to the 
string at the position of the chosen dislocation (e.g., along 
the arrow in Fig. [2]) . To relax the underlying translational 
order, the dislocation should move a distance of order re -1 
(which corresponds to the typical distance between grain 
boundary scars). This involves an energ y ba rrier whose 
order of magnitude estimated from Eq. (131 is YukT 1 . 



The final estimate for the structural relaxation time is 
thus 



1 



ZV 2 



exp 



E + E 1 (k<j)- 



(14) 



where D g is a "bare" glide diffusion constant (i.e. de- 
scribing small displacements of a free dislocation [12j). 
Eo is the energy barrier due to the Peierls potential and 
other subdominant terms, and E\ oc Ya 2 corresponds to 
the barrier estimated above. These predictions will be 
compared with simulation results below. 

The above arguments have been directly adapted from 
the theory of spherical crystals, with the necessary 
changes coming from the negative curvature. However, 
there may be a more serious difference between spheri- 
cal and hyperbolic geometries due to the fact that the 
atomic system is of finite extent in the former and possi- 
bly of infinite extent in the latter (with periodic bound- 
ary conditions, a proper thermodynamic limit can be 
taken). As mentioned above, bona fide thermodynamic 
phase transitions are possible in the hyperbolic plane 
H 2 , but not on the sphere. Strictly speaking, freezing 
of the 12 irreducible 5-fold disclinations on a sphere in 
an icosahedral superlattice only occurs at zero temper- 
ature. At finite temperature, provided of course that 
one waits long enough for equilibrium to be achieved, 
the 12 defects can move distances of the order of the 
sphere radius and the whole system continuously reor- 
ganizes. On H 2 on the other hand, one can envisage a 
'crystal of defects' formed by a periodic arrangement of 
the irreducible 7-fold disclinations (with or without as- 
sociated grain boundary scars) on a {3, 7} lattice. The 
latter lattice is the only one that is compatible with a 
predominance of 6-fold coordination numbers. Indeed, 
in the hyperbolic plane, each possible lattice possesses 
a uniquely determined lattice spacing (for a given cur- 
vature), which is characteristic of its symmetry. Thus, 
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there exists only one possible lattice compatible with 
the density of irreducible defects. Long range or quasi- 
long range order could then be achieved through a finite- 
temperature phase transition, much like the formation of 
an Abrikosov-like lattice of vortices in a uniformly frus- 
trated XY spin model in 2D flat space [4] . Checking this 
possibility in simulations of our atomic liquid model is 
however very hard for two reasons: the difficult equili- 
bration of the system at low temperature in the "glassy" 
liquid regime on one hand and the constraint due to the 
periodic boundary condition on the other. The latter 
constraint comes from the fact that the elementary cell 
of the periodic boundary condition procedure, with its 
12(g — 1) irreducible disclinations, should be compatible 
with the {3, 7} tiling of the hyperbolic plane with vertices 
centered on the irreducible defects [2~S1. 



IV. SIMULATION RESULTS AND DISCUSSION 

As already illustrated in Fig. 1, we find that as temper- 
ature is lowered, the number of defects decreases drasti- 
cally until one is left with an irreducible number of neg- 
ative 7-fold disclinations and an excess of dislocations 
formed by tightly bound pairs of 5-fold and 7-fold discli- 
nations. We show additional low-T atomic configurations 
in Fig. [3] for several values of the curvature and different 
periodic boundary conditions (octagonal cell with 12 ir- 
reducible negative disclinations and 14-gonal cell with 24 
such disclinations). Note that the scars centered on the 
excess disclinations only become evident for small kg. 

To provide a more quantitative test of the theoretical 
prediction, we have monitored as a function of temper- 
ature the number of 5-fold (i.e. q = +1) disclinations. 
Below a temperature roughly equal to the ordering tem- 
perature T* in flat space [25], these disclinations are all 
paired with 7-fold (i.e. q = —1) disclinations and, more- 
over, disclinations of higher charge magnitude \q\ > 2 
have almost completely disappeared. As seen in the plot 
shown in Fig. |1J the number of 5-7 dislocations decreases 
with decreasing temperature and seems to saturate at a 
nonzero value, which we extract through an exponential 
fit to the data. This saturation value, divided by the 
number of irreducible 7-fold disclinations in the cell (i.e. 
12(g — 1)), gives an estimate of the average number 
of dislocations which are contained at T = in a grain 
boundary scar radiating from an irreducible disclination. 
This allows us to compare the values obtained for several 
curvature magnitudes kg and periodic boundary condi- 
tions with the theoretical prediction in Eq. (11), only 



one adjustable parameter being needed for this exercise. 
The outcome is displayed in Fig. [5] and shows compatibil- 
ity between the simulation data and the prediction. For 
the smallest curvature parameter included, kg = 0.05, 
we are only able to provide an upper bound to Nd, as 
we cannot equilibrate the liquid at low enough tempera- 
ture with the available computer resources. On the other 
hand, the results for the two other curvatures are robust 



and the comparison between the octagonal and 14-gonal 
periodic boundary conditions show very little finite-size 
effects. 

From Fig. [5j one can see that the structure in grain 
boundary scars no longer remains stable above a thresh- 
old curvature magnitude k c g > 0.2. This threshold value 
is nonuniversal and depends on the type of interaction 
potential between the particles. (In Ref. [10] a simi- 
lar threshold value (a/R c ) ~ 0.2 has been found in an 
experimental system of colloidal particles on a sphere, 
but this may well be coincidental.) In addition, we 
have checked that the typical length of the arms ema- 
nating from an irreducible disclination is always of the 
order of, but smaller than, the theoretical cutoff radius 
r c = k^ 1 cosh - (7/6). Note that in the presently acces- 
sible temperature range, the system always remains in a 
liquid, i.e. disordered, phase, albeit very slowly relaxing 
and apparently on the verge of glass formation. We have 
never been able to observe crystallization of the lattice 
of defects discussed above. 

We have also investigated the low-T dy- 
namics of the system by monitoring the "self 
intermediate scattering function", F s (k,t) = 

( 1 /^ v )E^i( p -i+i^( cosh (^W)>: where dj(t) is 
the displacement of the jth atom during an elapsed 
time t, P_i +i h is a Legendre function of the first kind 
(which generalizes the standard function describing 
the correlation of density Fourier modes in 2D Eu- 
clidean space), and the wave- vector k is chosen as 
usual close to 2it/g [T7] [TB]. From F s (k,t), we have 
extracted the structural (translational) relaxation time 
r [T71 ITS] , which can then be studied as a function of 
temperature and curvature. One difficulty in trying to 
assess the validity of Eq. ( 14 ) is that the prediction 



only applies at low enough temperature in a regime 
dominated by irreducible grain boundary scars, when 
most thermal defects have disappeared and growth of 
the hexatic/hexagonal order has saturated. This occurs 
below a crossover temperature Tx, which in principle 
may be curvature-dependent. For the limited range over 
which we have simulation data (mostly, kg = 0.1 and 
0.2), Tx appears to vary weakly with curvature and we 
replace Eq. ( 14 ) by the following expression for T < Tx , 



r~r x exp( [E + E x (kg) l ] 



1 



1 



k B T k B T x 



, (15) 



where tx = t(Tx) depends on curvature. The latter de- 
pendence can be tentatively estimated by assuming that 
the relaxation at Tx is governed by the bare diffusion of 
dislocations to a distance of order k" 1 , i.e. tx — D^ 1 kT 2 . 

In Fig. [6j we test Eq. |l5} by plotting ln(r£» ff K 2 )/[l + 
B(kg) -1 ] versus 1/T for Toelow a Tx estimated around 
0.5T* and with the constants B and D g adjusted to pro- 
vide the best straight line. Notwithstanding several ad- 
justable parameters and a limited domain of available 
data, the dependence of the structural relaxation is at 
least compatible with the above prediction. However, a 
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FIG. 3: (Color online) Low-T atomic configurations showing the defect structures on H 2 : (a) k<t = 0.2, 
T/T* ~ 0.345, octagonal pbc; (b) k<j = 0.1, T/T* ~ 0.497, octagonal pbc; (c) kg = 0.1, T/T* ~ 0.520, 14-gonal pbc; 
(d) kit = 0.05, T/T* ~ 0.521, octagonal pbc. The system is shown at constant atomic size er, which allows one to 
better visualize the extension of the scars as a function of curvature. Note also the difference between the two 
different pbe's in (b) and (c): 12 irreducible defects are present for the octagon and 24 for the 14-gon. 
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FIG. 4: (Color online) Temperature dependence of the 
number of 5-fold disclinations divided by the number of 
irreducible 7-fold disclinations (i.e., 12 (g — 1)) for 
T < T*: no = 0.2, no = 0.1, no = 0.05. For no = 0.1 
and no = 0.2, both octagonal and 14-gonal pbc's are 
shown. The lines correspond to an exponential fit. 



more extensive data set would clearly be necessary to 
provide a more rigorous test of this finding. 



CONCLUSION 



The existence and the properties of these scars are in 
agreement with the theoretical predictions, which so far 
had only been confirmed for positively curved (spherical) 
geometry. Questions remain open for future work, such 
as an improved description of the long-time relaxation 
associated with large-distance diffusion of the defects and 
the possible occurrence of crystals of defects and scars. 
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In this article we have studied the low-temperature 
structure and dynamics of a Lennard- Jones liquid model 
embedded in manifolds of constant negative curvature 
and we have compared computer simulation results with 
predictions of a continuum theory of defects in curved 
space. In the regime where the spatial extension of the 
local hexagonal order has saturated due to the frustration 
effect, we find that the atomic configurations consist of 
isolated defect structures formed by an irreducible finite 
density of curvature-induced, negative disclination from 
which emerge finite-length strings of dislocations. Such 
"grain boundary scars" are unique to curved geometries. 



FIG. 5: (Color online) Extrapolated number of 
dislocations per irreducible disclination as a function of 

(kct) -1 . The straight line is the theoretical prediction 
from Eq. ( |Tlj ) (up to an adjustable parameter which is 
the intercept with the rr-axis). The data for kct = 0.05 

(octagonal pbc) is shown as an empty circle and only 
corresponds to an upper bound, as the liquid cannot be 
equilibrated at low enough T (see Fig. El) . 
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